CITEViz: interactively classify cell populations in CITE-Seq via a flow cytometry-like gating workflow using R-Shiny.

BACKGROUND: The rapid advancement of new genomic sequencing technology has enabled the development of multi-omic single-cell sequencing assays. These assays profile multiple modalities in the same cell and can often yield new insights not revealed with a single modality. For example, Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-Seq) simultaneously profiles the RNA transcriptome and the surface protein expression. The surface protein markers in CITE-Seq can be used to identify cell populations similar to the iterative filtration process in flow cytometry, also called "gating", and is an essential step for downstream analyses and data interpretation. While several packages allow users to interactively gate cells, they often do not process multi-omic sequencing datasets and may require writing redundant code to specify gate boundaries. To streamline the gating process, we developed CITEViz which allows users to interactively gate cells in Seurat-processed CITE-Seq data. CITEViz can also visualize basic quality control (QC) metrics allowing for a rapid and holistic evaluation of CITE-Seq data. RESULTS: We applied CITEViz to a peripheral blood mononuclear cell CITE-Seq dataset and gated for several major blood cell populations (CD14 monocytes, CD4 T cells, CD8 T cells, NK cells, B cells, and platelets) using canonical surface protein markers. The visualization features of CITEViz were used to investigate cellular heterogeneity in CD14 and CD16-expressing monocytes and to detect differential numbers of detected antibodies per patient donor. These results highlight the utility of CITEViz to enable the robust classification of single cell populations. CONCLUSIONS: CITEViz is an R-Shiny app that standardizes the gating workflow in CITE-Seq data for efficient classification of cell populations. Its secondary function is to generate basic feature plots and QC figures specific to multi-omic data. The user interface and internal workflow of CITEViz uniquely work together to produce an organized workflow and sensible data structures for easy data retrieval. This package leverages the strengths of biologists and computational scientists to assess and analyze multi-omic single-cell datasets. In conclusion, CITEViz streamlines the flow cytometry gating workflow in CITE-Seq data to help facilitate novel hypothesis generation.

Mutational screens highlight glycosylation as a modulator of colony-stimulating factor 3 receptor (CSF3R) activity.

The colony-stimulating factor 3 receptor (CSF3R) controls the growth of neutrophils, the most abundant type of white blood cell. In healthy neutrophils, signaling is dependent on CSF3R binding to its ligand, CSF3. A single amino acid mutation in CSF3R, T618I, instead allows for constitutive, ligand-independent cell growth and leads to a rare type of cancer called chronic neutrophilic leukemia. However, the disease mechanism is not well understood. Here, we investigated why this threonine to isoleucine substitution is the predominant mutation in chronic neutrophilic leukemia and how it leads to uncontrolled neutrophil growth. Using protein domain mapping, we demonstrated that the single CSF3R domain containing residue 618 is sufficient for ligand-independent activity. We then applied an unbiased mutational screening strategy focused on this domain and found that activating mutations are enriched at sites normally occupied by asparagine, threonine, and serine residues-the three amino acids which are commonly glycosylated. We confirmed glycosylation at multiple CSF3R residues by mass spectrometry, including the presence of GalNAc and Gal-GalNAc glycans at WT threonine 618. Using the same approach applied to other cell surface receptors, we identified an activating mutation, S489F, in the interleukin-31 receptor alpha chain. Combined, these results suggest a role for glycosylated hotspot residues in regulating receptor signaling, mutation of which can lead to ligand-independent, uncontrolled activity and human disease.

Disruption of the MYC Superenhancer Complex by Dual Targeting of FLT3 and LSD1 in Acute Myeloid Leukemia.

Mutations in Fms-like tyrosine kinase 3 (FLT3) are common drivers in acute myeloid leukemia (AML) yet FLT3 inhibitors only provide modest clinical benefit. Prior work has shown that inhibitors of lysine-specific demethylase 1 (LSD1) enhance kinase inhibitor activity in AML. Here we show that combined LSD1 and FLT3 inhibition induces synergistic cell death in FLT3-mutant AML. Multi-omic profiling revealed that the drug combination disrupts STAT5, LSD1, and GFI1 binding at the MYC blood superenhancer, suppressing superenhancer accessibility as well as MYC expression and activity. The drug combination simultaneously results in the accumulation of repressive H3K9me1 methylation, an LSD1 substrate, at MYC target genes. We validated these findings in 72 primary AML samples with the nearly every sample demonstrating synergistic responses to the drug combination. Collectively, these studies reveal how epigenetic therapies augment the activity of kinase inhibitors in FLT3-ITD (internal tandem duplication) AML. IMPLICATIONS: This work establishes the synergistic efficacy of combined FLT3 and LSD1 inhibition in FLT3-ITD AML by disrupting STAT5 and GFI1 binding at the MYC blood-specific superenhancer complex.

Asxl1 deletion disrupts MYC and RNA polymerase II function in granulocyte progenitors.

Mutations in the gene Additional Sex-Combs Like 1 (ASXL1) are recurrent in myeloid malignancies as well as the pre-malignant condition clonal hematopoiesis, where they are universally associated with poor prognosis. However, the role of ASXL1 in myeloid lineage maturation is incompletely described. To define the role of ASXL1 in myelopoiesis, we employed single cell RNA sequencing and a murine model of hematopoietic-specific Asxl1 deletion. In granulocyte progenitors, Asxl1 deletion leads to hyperactivation of MYC and a quantitative decrease in neutrophil production. This loss of granulocyte production was not accompanied by significant changes in the landscape of covalent histone modifications. However, Asxl1 deletion results in a decrease in RNAPII promoter-proximal pausing in granulocyte progenitors, indicative of a global increase in productive transcription. These results suggest that ASXL1 inhibits productive transcription in granulocyte progenitors, identifying a new role for this epigenetic regulator in myeloid development.

Integrative analysis of drug response and clinical outcome in acute myeloid leukemia.

Acute myeloid leukemia (AML) is a cancer of myeloid-lineage cells with limited therapeutic options. We previously combined ex vivo drug sensitivity with genomic, transcriptomic, and clinical annotations for a large cohort of AML patients, which facilitated discovery of functional genomic correlates. Here, we present a dataset that has been harmonized with our initial report to yield a cumulative cohort of 805 patients (942 specimens). We show strong cross-cohort concordance and identify features of drug response. Further, deconvoluting transcriptomic data shows that drug sensitivity is governed broadly by AML cell differentiation state, sometimes conditionally affecting other correlates of response. Finally, modeling of clinical outcome reveals a single gene, PEAR1, to be among the strongest predictors of patient survival, especially for young patients. Collectively, this report expands a large functional genomic resource, offers avenues for mechanistic exploration and drug development, and reveals tools for predicting outcome in AML.

GoPeaks: histone modification peak calling for CUT&Tag.

Genome-wide mapping of histone modifications is critical to understanding transcriptional regulation. CUT&Tag is a new method for profiling histone modifications, offering improved sensitivity and decreased cost compared with ChIP-seq. Here, we present GoPeaks, a peak calling method specifically designed for histone modification CUT&Tag data. We compare the performance of GoPeaks against commonly used peak calling algorithms to detect histone modifications that display a range of peak profiles and are frequently used in epigenetic studies. We find that GoPeaks robustly detects genome-wide histone modifications and, notably, identifies a substantial number of H3K27ac peaks with improved sensitivity compared to other standard algorithms.

PU.1 and MYC transcriptional network defines synergistic drug responses to KIT and LSD1 inhibition in acute myeloid leukemia.

Responses to kinase-inhibitor therapy in AML are frequently short-lived due to the rapid development of resistance, limiting the clinical efficacy. Combination therapy may improve initial therapeutic responses by targeting pathways used by leukemia cells to escape monotherapy. Here we report that combined inhibition of KIT and lysine-specific demethylase 1 (LSD1) produces synergistic cell death in KIT-mutant AML cell lines and primary patient samples. This drug combination evicts both MYC and PU.1 from chromatin driving cell cycle exit. Using a live cell biosensor for AKT activity, we identify early adaptive changes in kinase signaling following KIT inhibition that are reversed with the addition of LSD1 inhibitor via modulation of the GSK3a/b axis. Multi-omic analyses, including scRNA-seq, ATAC-seq and CUT&Tag, confirm these mechanisms in primary KIT-mutant AML. Collectively, this work provides rational for a clinical trial to assess the efficacy of KIT and LSD1 inhibition in patients with KIT-mutant AML.

Mutated SETBP1 activates transcription of Myc programs to accelerate CSF3R-driven myeloproliferative neoplasms.

Colony stimulating factor 3 receptor (CSF3R) mutations lead to JAK pathway activation and are the molecular hallmark of chronic neutrophilic leukemia (CNL). Approximately half of patients with CNL also have mutations in SET binding protein 1 (SETBP1). In this study, we developed models of SETBP1-mutated leukemia to understand the role that SETBP1 plays in CNL. SETBP1 mutations promote self-renewal of CSF3R-mutated hematopoietic progenitors in vitro and prevent cells from undergoing terminal differentiation. In vivo, SETBP1 mutations accelerate leukemia progression, leading to the rapid development of hepatosplenomegaly and granulocytosis. Through transcriptomic and epigenomic profiling, we found that SETBP1 enhances progenitor-associated programs, most strongly upregulating Myc and Myc target genes. This upregulation of Myc can be reversed by LSD1 inhibitors. In summary, we found that SETBP1 mutations promote aggressive hematopoietic cell expansion when expressed with mutated CSF3R through the upregulation of Myc-associated gene expression programs.

Co-Occurring CSF3R W791* Germline and Somatic T618I Driver Mutations Induce Early CNL and Clonal Progression to Mixed Phenotype Acute Leukemia.

Chronic neutrophilic leukemia (CNL) relates to mutational CSF3R activation with membrane proximal CSF3R mutations such as T618I as driver mutations, but the significance of truncating mutations is not clarified. In CNL, concomitant mutations promote disease progression, but insight into longitudinal acquisition is incomplete. In this study, we investigated the role of co-occurring germline and somatic CSF3R mutations in CNL, and assessed the impact of clonal evolution on transformation to acute leukemia. We employed sequential next generation sequencing and SNP array karyotyping to assess clonal evolution in CNL of early manifestation age based on a 33-year-old patient. Germline vs. somatic mutations were differentiated using a sample from the hair follicle. To investigate a potential predisposition for CNL development and progression by germline CSF3R-W791*, allelic localizations were evaluated. We detected a somatic CSF3R-T618I mutation at 46% variant allele frequency (VAF) at the time of CNL diagnosis, which co-occurred with a CSF3R-W791* truncation at 50% VAF in the germline. Evaluation of allelic localization revealed CSF3R-T618I and W791* on the same allele. A concomitant ASXL1 mutation at 39% VAF increased to 48% VAF upon transformation to mixed phenotype acute leukemia (MPAL), which has both myeloid and lymphoid features. Clonal evolution further involved expansion of the CSF3R double-mutant clone to 90% VAF via copy neutral loss of heterozygosity on chromosome 1p and the emergence of a RUNX1 mutant subclone. Allogeneic transplantation induced complete remission. This study highlights that CNL not only transforms to AML but also to MPAL. The molecular evolution is especially interesting with a CSF3R-W791* mutation in the germline and acquisition of CSF3R-T618I on the same allele compatible with increased susceptibility for mutation acquisition facilitating RUNX1-related clonal transformation.

TNK1 is a ubiquitin-binding and 14-3-3-regulated kinase that can be targeted to block tumor growth.

TNK1 is a non-receptor tyrosine kinase with poorly understood biological function and regulation. Here, we identify TNK1 dependencies in primary human cancers. We also discover a MARK-mediated phosphorylation on TNK1 at S502 that promotes an interaction between TNK1 and 14-3-3, which sequesters TNK1 and inhibits its kinase activity. Conversely, the release of TNK1 from 14-3-3 allows TNK1 to cluster in ubiquitin-rich puncta and become active. Active TNK1 induces growth factor-independent proliferation of lymphoid cells in cell culture and mouse models. One unusual feature of TNK1 is a ubiquitin-association domain (UBA) on its C-terminus. Here, we characterize the TNK1 UBA, which has high affinity for poly-ubiquitin. Point mutations that disrupt ubiquitin binding inhibit TNK1 activity. These data suggest a mechanism in which TNK1 toggles between 14-3-3-bound (inactive) and ubiquitin-bound (active) states. Finally, we identify a TNK1 inhibitor, TP-5801, which shows nanomolar potency against TNK1-transformed cells and suppresses tumor growth in vivo.

Mutant SETBP1 enhances NRAS-driven MAPK pathway activation to promote aggressive leukemia.

Mutations in SET-binding protein 1 (SETBP1) are associated with poor outcomes in myeloid leukemias. In the Ras-driven leukemia, juvenile myelomonocytic leukemia, SETBP1 mutations are enriched in relapsed disease. While some mechanisms for SETBP1-driven oncogenesis have been established, it remains unclear how SETBP1 specifically modulates the biology of Ras-driven leukemias. In this study, we found that when co-expressed with Ras pathway mutations, SETBP1 promoted oncogenic transformation of murine bone marrow in vitro and aggressive myeloid leukemia in vivo. We demonstrate that SETBP1 enhances the NRAS gene expression signature, driving upregulation of mitogen-activated protein kinase (MAPK) signaling and downregulation of differentiation pathways. SETBP1 also enhances NRAS-driven phosphorylation of MAPK proteins. Cells expressing NRAS and SETBP1 are sensitive to inhibitors of the MAPK pathway, and treatment with the MEK inhibitor trametinib conferred a survival benefit in a mouse model of NRAS/SETBP1-mutant disease. Our data demonstrate that despite driving enhanced MAPK signaling, SETBP1-mutant cells remain susceptible to trametinib in vitro and in vivo, providing encouraging preclinical data for the use of trametinib in SETBP1-mutant disease.

Outgrowth of a CSF3R-mutant clone drives a second myeloproliferative neoplasm in a chronic myeloid leukemia patient: a case report.

BACKGROUND: Chronic myeloid leukemia (CML) and chronic neutrophilic leukemia (CNL) are two myeloproliferative neoplasms with mutually exclusive diagnostic criteria. A hallmark of CML is the Philadelphia chromosome (Ph), which results in a BCR-ABL1 fusion gene and constitutive tyrosine kinase activity. CNL is a Ph-negative neoplasm and is defined in part by the presence of CSF3R mutations, which drive constative JAK/STAT signaling. CASE PRESENTATION: Here, we report the exceedingly rare co-occurrence of two granulocytic myeloproliferative neoplasms in a 69-year old male patient. After an initial diagnosis of chronic myeloid leukemia, the patient's clinical course was shaped by hematologic toxicity, the emergence of treatment-resistant BCR-ABL1 clones, and the expansion of a CSF3R-mutant clone without ABL1 mutations under selective pressure from tyrosine kinase inhibitors. The emergence of the CSF3R-mutant, neutrophilic clone led to the diagnosis of CNL as a second myeloproliferative neoplasm in the same patient. CONCLUSIONS: This is the first reported case of CNL arising subsequent to CML, which occurred under selective pressure from targeted therapy in a patient with complex clonal architecture. Patients with such molecularly complex disease may ultimately benefit from combination therapy that targets multiple oncogenic pathways.

Combined inhibition of JAK/STAT pathway and lysine-specific demethylase 1 as a therapeutic strategy in CSF3R/CEBPA mutant acute myeloid leukemia.

Acute myeloid leukemia (AML) is a deadly hematologic malignancy with poor prognosis, particularly in the elderly. Even among individuals with favorable-risk disease, approximately half will relapse with conventional therapy. In this clinical circumstance, the determinants of relapse are unclear, and there are no therapeutic interventions that can prevent recurrent disease. Mutations in the transcription factor CEBPA are associated with favorable risk in AML. However, mutations in the growth factor receptor CSF3R are commonly co-occurrent in CEBPA mutant AML and are associated with an increased risk of relapse. To develop therapeutic strategies for this disease subset, we performed medium-throughput drug screening on CEBPA/CSF3R mutant leukemia cells and identified sensitivity to inhibitors of lysine-specific demethylase 1 (LSD1). Treatment of CSF3R/CEBPA mutant leukemia cells with LSD1 inhibitors reactivates differentiation-associated enhancers driving immunophenotypic and morphologic differentiation. LSD1 inhibition is ineffective as monotherapy but demonstrates synergy with inhibitors of JAK/STAT signaling, doubling median survival in vivo. These results demonstrate that combined inhibition of JAK/STAT signaling and LSD1 is a promising therapeutic strategy for CEBPA/CSF3R mutant AML.